﻿ 6 Analiza Fourier a semnalelor definite in timp discret 1-Seria Fourier pentru semnale periodice in timp discret 2-Transformata Fourier in timp discret pentru semnale aperiodice http://shannon etc upt ro/teaching/ssist/Cap6 pdf 1 Raspunsul sistemelor discrete, liniare si invariante in timp SLITD la exponentiala complexa de modul unitar 0jn jn  0eH() 0eH  Functie proprie pentru Valoare proprie jk H h k e        orice SLITD k   •Demonstratie:   kjnjknjnj  0000ekheekhenhny  kk         00jnjnjk  0                     notatia: H h k e y n e H H e00Cu k   •La iesire se obtine tot o exponentiala complexa, de aceeasi frecventa ca si semnalul de intrare, , dar de modul si faza 0 initiala afectate de |H()| si arg{H()} 00 2 1 Fie un semnal de intrare x[n] ce apare ca o combinatie liniara a unor exponentiale complexe de forma Φ[n]=exp(jn), kk Se observa ca raspunsul este o combinatie liniara a raspunsurilor partiale, H(Ω) Φ[n] kk kjn aekjn    kka H e kH() kk  Daca semnalul de intrare se poate pune sub forma unei sume de exponentiale discrete, este suficienta cunoasterea functiei H(), pentru a determina raspunsul y[n] 3 Seria Fourier in timp discret pentru semnale discrete si periodice Fie un semnal periodic, x[n] de perioada N, pentru care x[n+N]=x[n], n Intr-o perioada exista N valori {x , x , …, x[N-1]}, dupa care ele se repeta: x = x[N], x[N+1]= x etc Se scrie ca: x[n]= x[(n)] N (n)- reprezentarea lui n in clase de resturi modulo N N (7)=7 pentru ca 7=08+7 8 (15)=7 pentru ca 15=18+7 8 Pentru n •Periodicitatea seriei Fourier se verifica usor: 2211NN j k N n jk n    211jn NN        kNc x n e x n e e 00nnNN 21N jk n 21jn N Dar 1 si , 0 1 e c x n e c k N         k N k 0nN c c c c   sau k N k kk   N 6 3 Exemple 2      nnxsin   1 Semnalul sinusoidal are perioada N Folosind N   relatia lui Euler:    2222      111121nNjnjnjnj      sinNNNNeeeen   2222jjjjN  Prin identificare: 11        1 1 0 2 2 , ; in rest coeficientii sunt nuli: 0 NNc c c c c 22jj sin 2 / 6n  Diagrame spectrale de modul si de faze pentru N=6  7 2 1 sin 2 / 4cos 2 / cos 4 / / 2x n n N n N n N                 Relatia lui Euler:  12211Nj n j n  1NNx n e e    22jj  122Nj n j n  2 2NNee   222211Nj n j njj   22NN e e e e 22 Identificare: 11jj         0 1 2 2 11; 2 ; ; 2NNc c c c c 2 2 2 2jj 8 4 3 Semnalul dreptunghiular, de perioada N Intr-o perioada sunt 2N+1 valori succesive egale cu 1 1 n   1122 2 2NN  111jk n jk N jk ; 0 1N N Nc e e e k N      k  NN   10n N n 21N Pentru k=0: 1c 0 N 1 1kN       1212N   1sinsin 2 1kN  1122N    c     2kNN  sinsink 22N   k N 9 •Semnale in timp continuu: trunchierea seriei duce la aparitia unor oscilatii la trecerea brusca a semnalului de la o valoare la alta Mariind numarul de termeni din serie, amplitudinea oscilatiilor nu se modifica, dar ele devin mai rapide (fenomenul Gibbs) •Semnale in timp discret: trunchierea nu duce la aparitia unui fenomen de tip Gibbs Se face o aproximare a semnalului; cu cat numarul de termeni e mai mare, cu atat aceasta este mai buna Cand se insumeaza toti cei N termeni, semnalul este chiar x[n], fara nici o eroare •Fie N = 9 si 2N+1 = 5 Coeficientii Fourier sunt: 1 c = 0,556 , 0 c = c =0,32 , 18 c = c = -0,059 , 27 c = c = -0,111, 36 10 c = c = 0,073 45 5 Semnalul trunchiat cu 2M+1 termeni din seria sa Fourier: 2M jk n 9    Mkx n c e nM Pentru M=1, 2, 3 si 4 avem aproximarile: 2n   10 556 0 64 cosxn 9 24nn     20 556 0 64 cos 0 118cosxn 99 246nnn      30 556 0 64 cos 0 118cos 0 222 cosxn 9 9 9 2 4 6 8n n n n          40 556 0 64 cos 0 118cos 0 222 cos 0 146 cosxn 9 9 9 9 11 Trunchierea seriei de reconstructie pentru semnalul dreptunghiular Nu apare o comportare asemanatoare fenomenului Gibbs: 12 x[n]=x[n] 4 6 Proprietatile seriei Fourier in timp discret 1 Liniaritatea x[n], y[n] periodice de aceeasi perioada yx           kkax n by n ac bc Tema: demonstratia 2 Deplasarea (translatarea) in timp 2  0jk n N      0kx n n e c    Spectrul de modul nu este afectat, in schimb este afectat cel de faza Operatii duale: deplasarea si modulatia (inmultirea cu o exponentiala complexa) 13 3 Conjugarea complexa  ** x n c c        kk   N   4 Reflectarea semnalului            kkx n c c    N Reflectarea este auto-duala 14 7 5 Modificarea scarii timpului Semnal cu derulare de m ori mai lenta: / ; daca divisibil cu 1x n m n m n m        , perioada 'x n c N mN       mk    0 ; altfel m  Semnalul x[n] periodic de perioada N => Se obtine semnalul x[n] periodic cu perioada N’=mN (m) 2211mN N  11jk jk mpn mN mN'      mmkc x n e x pm e     00npmN mN 21N  11jk p Nx p e c   k 0pmN m 15 6 Modularea semnalului Modularea realizeaza deplasarea cu k a spectrului de modul si faza 0 2 0jkn N e x n c c         0kkkk   0 N 7 Produsul a doua semnale (teorema produsului) Pentru x[n], y[n] periodice de aceeasi perioada 1N yx x x         mkkx n y nc c c c km    0Nm  Seria Fourier este convolutia periodica sau circulara a secventelor discrete formate din cx si cy kk 16 8 8 Convolutia periodica (teorema convolutiei) x[n] si y[n] periodice de perioada N, cu seriile Fourier cx si cy kk Convolutia lor circulara este z[n], de asemenea periodica cu N: 1N    z n z n N                Nz n x n y n x k y n k  0k 1N x n y n x m y n m          0m 1N x m y n m      N  0m yx x n y n Ncc         kk 17 9 Diferentierea discreta a semnalului discret Echivalenta diferentierii in timp continuu 2   jk  11Nx n x n e c           k     Tema: demonstratie 10 Insumarea in domeniul (timpului) n Suma esantioanelor unui semnal x [n] de perioada N , fara componenta continua ( ) 0care aceeasi perioada N 0  nc k       02, 0y n x m c jk   mN 1e  Demonstratie 2  jk  yxN 11y n y n x n e ccc           k k k 18   9 11 Proprietati specifice semnalelor reale Daca x[n] este real,            kkkx n x n c c c    N   ; Arg Arg Argc c c c c c     k k k kkk      NN  Re Re ; Im Im Imc c c c c            k k k kk    N 12 Seria Fourier a componentei pare si impare     11x n x n       Rex n c c c        pk k k 2 2 2     11x n x n       Imx n c c j c        ik k k 2 2 2 19 13 Relatia lui Parseval 11NN 222 x n x n N c      2k 00nk  Patratul normelor in l2 Puterea P a semnalului discret periodic calculabila -in timp: ca medie a patratelor esantioanelor temporale, -in frecventa: ca suma a patratelor modulelor coeficientilor seriei Fourier atasate semnalului 11NN 221   kP x n c 00nkN 20 10 Transformarea Fourier in timp discret, pentru semnale discrete Introducem notiunea de transformata Fourier, pornind de la seria Fourier a unui semnal periodic, a carui perioada tinde la infinit N    x n x n kN      k   Semnalul de durata finita este: Pentru semnalul periodic avem: 2 2 jk n1jk n NN c x n e      kx n c e k kN   kNN   Seria Fourier se poate rescrie: 2  1jk n N c x n e  k21 nN Se defineste anvelopa produsului Nck prin:  jn X x n e      n   Coeficientii semnalului periodizat devin : 12  c X k,     00k NN Pentru semnalul considerat    1sin 2 1N 2  X   sin 222 11 2N+1=5 1 Spectrul semnalului neperiodic x[n] este anvelopa sau infasuratoarea X(Ω) 23 Relatia de recuperare a semnalului este: 11jk n jk n 00 x n X k e X k e           0 0 0 2N k   N  k   N  La limita N  (Ω 0) , obtinem semnalul neperiodic: 0 1 jn  limx n x n X e d           2N 2   jn X x n e   Transformata Fourier in timp discret  n 1jnTransformata Fourier inversa in timp x n X e d      2 2discret 12     kc X k00, NN24 12 Proprietati •Pentru semnale absolut sumabile, x[n]  l1 , transformata converge    jnjn                    nxnxenxenxX   1     nnn  •Functia X() este continua X(Ω) – spectrul unui semnal discret este continuu •Spectrul unui semnal discret este o functie periodica dupa , cu perioada 2 Este deci suficienta reprezentarea pe intervalul [0,2π) sau [-π, π)    2j n j n   2X x n e e X            n   2XX                 22   njnjnjnj    eeee 25 2 x n l Semnale de energie finita (semnale de patrat sumabil)   Convergenta in medie patratica N jn l i mX x n e      N nN N jn lim 0X x n e        N 2nN Trunchierea seriei duce la aparitia fenomenului Gibbs (ca si in timp continuu) Probleme de convergenta apar doar la semnalele cu suport infinit 1 Semnalele cu suport finit sunt din l 26 13  Exemple  s in 2 N 1   1 2       1x n n N n N X                 1 11  sin 2 2 Semnalul treapta unitate  0jn          1x n n X n e e          n   1n         3 Semnalul cauzal n  , 1x n a n a      Functie complexa, modul si faza 1 sina  X arctg;        2 cos 1a 1 2 cosaa   27 1a sin   X Arg X arctg ;          21a cos 12a cos a    28 14 Transformarea Fourier in timp discret, pentru semnale discrete si periodice •Poate fi introdusa numai in sens distributional  2   X c l c c2 , unde          ll l N     NlN    Un semnal discret periodic, poate fi descompus in seria N 1  njk    ,0 kecnx k 0 Exponentiala complexa cu frecventa kΩ de modul unitar are 0 transformarea in timp discret: 0jk n          02ek 29 2 Folosind liniaritatea avem 1N             02kx n c k 2 0k Transformarea exponentialei complexe de modul unitar  0jn                     002 2 2n e X k 2 k 30 15 1N             02kx n c k 2 0k Tinand seama de definitia distributiei Dirac periodice: 11NN    2                               02 2 2 2kkX c k m c k m          00mmkkN   Fie l = k+mN Tinand cont de periodicitatea seriei Fourier putem scrie c= c l k  2   X c l c c2 , where         ll l N      NlN    Am obtinut transformarea Fourier in timp discret pentru un semnal 31 periodic discret Exemplu: Distributia Dirac periodica  1 n n kN c              seria FourierNk kN    2  2X c l         l    lN    2 2 2    nk         ;     000N  kN N N 32 16 Spectrul X() unui semnal discret periodic este format din   distributii Dirac plasate la 2/kN k Spectrul este periodic de perioada 2 33 Proprietati 1 Liniaritatea ax n by n aX bY              2 Deplasarea in timp a semnalului 0jn       0x n n e X  3 Modularea in domeniul timp 0jn      0e x n X  4 Scalarea variabilei timp x n X k      k   34 17 5 Conjugarea complexa a semnalului ** x n X        6 Reflectarea in timp a semnalului x n X         7 Diferentierea numerica a semnalului discret j  11x n x n e X             8 Convolutia semnalelor (teorema convolutiei) x n y n X Y             35 9 Insumarea in domeniul timpului n X               20x k X  j 1e k    10 Produsul semnalelor discrete 11    x n y n X u Y u du X Y               22     2 11 Derivarea in domeniul spectrului dX  nx n j   d 36 18 12 Spectrele semnalelor reale ** x n x n X X               x n Re X x n j Im X  ;               pi     X X Arg X Arg X ; ;                     ; Im XRe X Re X Im X                       13 Relatia lui Parseval 221 2 ;x n l , x n X d        2 2n- 22     222X x n  ,Ll    2        X 2x n , y n l , ,Y x n , y n22 Ll37 ,    Densitatea spectrala de energie Se noteaza cu S(Ω) X 2  S X ;     X  1     XW S d   2 2 Energia se obtine prin integrarea densitatii spectrale de energie 38 19 Raspunsul in frecventa al SLITD Raspunsul la impulsul unitar: h[n] h n H      Transformata sa Fourier este raspunsul sistemului in frecventa H() Cunoscand H() se poate afla iesirea pentru orice intrare i)Se determina spectrul semnalului de intrare x[n], X(), ii)Se determina spectrul semnalului de iesire y[n], Y(), iii) Se descompune Y() intr-o suma de fractii simple si se aplica transformarea Fourier inversa Pentru aceasta, se 39 folosesc de obicei tabelele de transformate! Raspunsul unui SLITD la un semnal de intrare discret si periodic 11NN 2jk n jk n  00             kkx n c e y n c H k e 00 ; 00kkN In cazul particular semnal de intrare armonic: 2                0 0 0 00cos cosx n A n y n A H n   N 40 20 Demonstratie!    22          00jjA     2NN                0cosx n A n e e 2N   Coeficientii Fourier:   00jjAA c e c c e;   1 1 1N 22       0 0 0 02nnj j j jAA             y n e H e e H e0 0 0; 22N     H H h n fiindca                    0 0 00 0 0jnjnAA            00y n H e H e 22 2 y n A H n         cos ;    0 0 0 0 041 N SLITD descrise de ecuatii cu diferente finite liniare cu coeficienti constanti Forma generala a ecuatiei este NM         0, 0kka y n k b x n k a 00kk Se aplica transformarea Fourier : Rezulta raspunsul in frecventa al sistemului 42 21 Exemple i) 21  1 2 1y n y n y n x n x n                24 Raspunsul in frecventa jj    11ee H     jj         22111jjee      12  1ee 24  j 421      1,21je 42 Se face descompunerea in fractii simple: AA 121 2 2 1   H A j;     1,2jj    2211ee    12 43 Luand transformarea inversa se obtine raspunsul la impuls al sistemului nn           1 1 2 2h n A A n   jn jn     441 1 2 2 1 1 2 2 1 j e j e n           n  2 2 2 2 2       1  cos 2 2 1 sinn n n       n   2 4 4 21n  cos 2sinnn      n  2 4 4 44 22 ii) Determinam raspunsul la impuls al sistemului din raspunsul in frecventa: 3 H   11jj        11ee     28    Descompunere in fractii simple: 41 H     11jj    11ee 28 Transformarea Fourier inversa : 11 4h n n       3nn 22 45 iii) Obtinem raspunsul la impuls h[n] si raspunsul la treapta unitara (functia indiciala) s[n] pentru sistemul descris de ecuatia  1 , 1y n ay n x n a         Se identifica coeficientii si avem: 1 n   H h n a n raspunsul la impuls         j 1ae Transformarea Fourier a semnalului treapta unitara σ[n] este: 1 X      02j  1e Raspunsul la treapta unitara, s[n], are transformarea: 1 S H X         02jj     111aae e   46 23 Descompunere in fractii simple: 1 1 1a S            2jj     1 1 111a a aae e    1 1 1a S            2jj      1111aaae e   Se aplica transformarea Fourier inversa si rezulta functia indiciala s[n] 1n 11aa n s n a n n n               1 1 1a a a   47 SLITD de ordinul intai si doi Raspunsul in frecventa a unui SLITD este o functie rationala Numitorul si numaratorul pot fi scrise ca si produs de factori cu coeficienti reali de gradul unu si doi 2P M P 2j j j                1211k k ke e e 0 1 1kkb H   2Q N Q   20j j ja              1211k k ke e e 11kk Implementarea : sisteme de ordinul unu sau doi conectate in serie realizate in forma directa II (Produsul in domeniul frecventa inseamna convolutie in 48 timp) 24 2P M P   2j j j              1211k k ke e e 0 1 1kkb H   2Q N Q   20j j ja              1211k k ke e e 11kk Pentru M=N, descompunerea in fractii simple  2jQ N Q   eAb 01N k k k H         2j j j          1111kkNae e e 12k k k Implementarea : sisteme de ordinul unu sau doi conectate in paralel realizate in forma directa II 49 •Un polinom de ordinul intai  sistem de ordin intai •Un polinom de ordinul doi  sistem de ordin doi Cunoscand comportamentul in frecventa al sistemelor de ordinul intai si doi, putem deduce comportamentul in frecventa oricarui sistem SLITD Rezulta importanta sistemelor de ordinul unu si doi 50 25 Sisteme discrete de ordinul intai 1  y n ay n x n a H1 , 1 ;             j 1ae 1n 1a n  h n a n s n n            1a Pentru 0<a<1 raspunsul la impuls precum si raspunsul indicial nu prezinta oscilatii 51 Pentru -1<a<0 raspunsurile au caracter oscilant! 1  y n ay n x n a H1 , 1 ;             j 1ae 1n 1a n  h n a n s n n            1a 52 26 Sisteme discrete de ordinul doi Ecuatia cu doi parametri: Raspunsul in frecventa: Raspunsul la impuls: Raspunsul indicial: 53 Sistem de ordinul doi: =0 1 0  1n 1 H h n n n         2 j 1e  1 1 1 1            2S  2 2 2 2jj    111jee  1 1 1       1e   1nn 1s n n n          22 1 11      54 27 =  2 Sistem de ordinul doi:     1 H    2 j 1e   n  1 ,h n n n             1nn       1s n n n               22 1 11          55 Sistem de ordin doi Raspunsul la impuls pentru =/2, =0, = Sistem critic amortizat Sistemul prezinta oscilatii 56 28 Raspunsul in frecventa 1 H    2 2 2 2 2 2 2 4 4 1 1 4cos cos cos cos                  2r sin a r cos        arctg    2 1 2 2r cos cos r cos     57 Functia de corelatie Densitatea spectrala de putere si de energie a semnalelor discrete N 1   xR k lim x* n x n k  21NN nN x[n] este de perioada N:periodic0Daca 01N 2 jm k 02     xmR k c e , 0 00nN  xR k N este periodica de perioada 0Functia 0011NN 02jm k          x x mR k R k N x* n x n k c e ,0 N 00nn Teorema Wiener-Hincin: are coeficientii seriei Fourier exponentialeR x 2 R k c  xm58 29 01N 2 R P c ; R k P0       x m x  0n - puterea medie pe o perioada 0P a semnalului x[n]N R k R k     xx Se defineste spectrala de putere a semnalului x[n]densitatea ca transformata Fourier in timp discret a semnalului periodic: 01N 2                 x x mR k S c k 02 0n o alta forma a teoremei Wiener-Hincin 59 Pentru un semnal de energie infinita, dar putere medie finita, neperiodic, densitatea spectrala de putere este transformata Fourier a functiei de corelatie:       xxS R k periodica, perioada 2 , nenegatS   xiva, para 1 0P R S d       xx   2 2 60 30 Functia de intercorelatie pentru semnale discrete de energie infinita dar de putere medie finita N 1        xyR k lim x* n y n k  21NN  nN cazul a doua semnale periodice de aceasi perioada N:0In     0xy xyR k N R k  Functia de intercorelatie a semnalelor discrete, periodice, aceeasi perioada: 01N 1        xyR k x* n y n k  00nN 61 Proprietati * R k R k   xy yx 2 00R k R R P P   xy x y x y puterile semnalelor discrete x[n], y[n]P ,P xy Pentru x y R R R   functia de intercorelatie xy xx x devine functia de corelatie sau autocorelatie 62 31 Functia de corelatie pentru semnale discrete de energie finita Masoara gradul de asemanare a celor doua semnale         xyR k x* n y n k n   x* k y k x * k y k                 - densitate interspectrala de energie,R k X * Y S         xy XY  x n y n , functia de autocorelatie                    xR k x* n x n k x k x k  n   63 Proprietati ale functiei de autocorelatie a semnalelor discrete de energie finita Teorema Wiener-Hincin 2          xxR k X S  21 0W R X d        xx   2 2 Functia este para si isi atinge maximul in origine Rk   x 64 32 65 Relatia intre densitatile spectrale de putere si de energie ale semnalelor ce trec prin sisteme discrete, linare si invariante in timp y n x n h n   2     YXS H S  XYS sau densitate spectrala de putere (pentru semnalele de putere medie finita) sau densitate spectrala de energie (pentru semnalele de energie finita)  y x hR n R n R n 66 corelatoare 33 